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The search for heavy Higgs bosons is an essential step in the exploration of the Higgs sector and 
in probing the Supersymmetric parameter space. This paper discusses the constraints on the Ma 
and tan /3 parameters derived from the bounds on the different decay channels of the neutral H and 
A bosons accessible at the LHC, in the framework of the phenomenological MSSM. The implications 
from the present LHC results and the expected sensitivity of the 14 TeV data are discussed in terms 
of the coverage of the [Ma — tan 0] plane. New channels becoming important at 13 and 14 TeV for 
low values of tan/3 are characterised in terms of their kinematics and the reconstruction strategies. 
The effect of QCD systematics, SUSY loop effects and decays into pairs of SUSY particles on these 
constraints are discussed in details. 

PACS numbers: 12.60. Jv, 14.80.Da 



I. INTRODUCTION 

With the observation of a light Higgs-like particle by 
the ATLAS and CMS experiments at the LHC [H |2], 
the detailed exploration of the Higgs sector becomes one 
of the most compelling programs of collider physics. In 
particular, understanding whether this sector extends be- 
yond that of the Standard Model (SM) and heavier Higgs 
bosons exist is of crucial importance for the viability 
of several models of new physics beyond the Standard 
Model, in primis of Supersymmetry (SUSY). This ques- 
tion can be answered either through a precision study 
of the couplings of the lightest boson, h, or by direct 
searches of the additional, heavier states which charac- 
terise extended Higgs models. 

The LHC experiments have not only observed a light 
state and obtained the first determination of its decay 
rates to 77, WW and ZZ. They have also performed 
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several searches directly probing the possible production 
of heavy Higgs bosons and other searches, which can now 
be re-interpreted in order to set constraints on the pro- 
duction and decays of neutral heavy Higgs states. How- 
ever, these data are still largely fragmentary. 

Several studies of the MSSM heavy Higgs sector at 
LHC results have already been performed [3HZ]. This 
paper intends to provide a comprehensive assessment of 
the present status and the future perspectives for the con- 
straints on the MSSM Higgs sector parameters, from the 
identification of the main processes relevant to the LHC 
searches to a systematic study of the exclusion limits de- 
rived from the combination of the LHC results, in the 
context of the phenomenological MSSM (pMSSM) with 
the neutralino as the lightest SUSY particle (LSP) [8]. 
We perform this study taking the mass of the heavy pseu- 
doscalar, Ma, and the ratio of the vacuum expectation 
value of the two Higgs doublets, tan /?, as the main pa- 
rameters. We highlight the complex pattern of decays 
arising at low values of tan/3, values which are shown 
to be compatible with the present data and discuss the 
complementarity of decay modes such as H — > ZZ, tt 
and hh and A — >• Zh. 

The combination of the relevant decay channels to ex- 
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tend the sensitivity of the heavy Higgs searches over most 
of the [Ma — tan/3] was already discussed in [9 . Here, 
we use the published and preliminary results for the ex- 
pected upper limits on the product of production cross 
section and decay branching fraction in several channels 
as constraints and extrapolate them to the full 2012 data 
set of 25 fb _1 /experiment at 8 TeV and to 150 fb _1 of 
14 TeV data. 

In section II, we discuss the production and decays 
of the H and A neutral bosons, with special empha- 
sis for the low tan /3 region. Section III is devoted to 
presentation of the results of our systematic study of 
the indirect constraints derived from the latest measure- 
ments of the decay rates of the 126 GeV Higgs-like parti- 
cle together with those obtained from direct searches for 
H/A ->> t + t", H S m ZZ, bbH ->> bbbb and resonant tt 
production and their expected sensitivity on the 14 TeV 
LHC data. Then, we review additional decay modes, 
which have not yet been considered in the LHC searches 
but will become important at 14 TeV in the low tan/3 
region and characterise their kinematics and reconstruc- 
tion strategies. Finally, we discuss the validity of these 
bounds when taking into account the production cross 
section uncertainties and the role of SUSY particles af- 
fecting the decays of heavy Higgs bosons, either in their 
direct decays to SUSY states or through loop corrections 
to their decay widths. Section IV has the conclusions. 



II. THE HIGGS SECTOR AND THE M A - tan/3 
PARAMETERS 

A. H and A production and decays in the pMSSM 

The MSSM neutral heavy Higgs bosons H and A have 
couplings modified compared to the SM Higgs state. In 
the decoupling limit (Ma ^ Mz), the H/A coupling to 
the top quarks is suppressed by 1/ tan/3, while the cou- 
plings to bottom quarks and tau leptons are enhanced by 
tan/3. As a consequence, the H/Att coupling is impor- 
tant only for tan/3 < 10, those to bb and tt becoming 
dominant for larger values. On the other hand, the H 
couplings to vector bosons are suppressed by a factor 
cos(/3 — a), which in the large Ma and tan/3 limit, de- 
creases as 1/ tan /3. The situation is the same for the AhZ 
coupling, while there is no A coupling to vector bosons 
at tree level. Finally, the coupling of the H to hh also 
decreases in the large Ma and tan/3 limit with 1/ tan/3. 
Hence, the description of the heavy Higgs sector in the 
large tan /3 limit is simply dominated by the couplings to 
b and r fermions, whereas in the small tan /3 regime, a 
rich phenomenology emerges, as the other couplings be- 
come important. A thorough discussion can be found in 
Ref. [10]. 

The H and A production cross section is dominated 
by the gluon fusion process and the associate Higgs pro- 
duction with b quarks. The relevant cross sections are 
shown in Fig. [l] for two values of the pseudoscalar Higgs 
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FIG. 1: Production cross sections for H and A bosons as a 
function of tan /3 in pp collisions at 8 TeV (right panels) and 
14 TeV (left panels) for M A = 300 GeV (upper panels) and 
500 GeV (lower panels). 



mass (M A = 300 and 500 GeV) at y/s = 8 and 14 TeV, 
as a function of tan /?. 

The bbH associate production is a tree level process, 
which increases as tan 2 /3 and becomes dominant for 
tan/3 > 10. Instead, the gluon fusion processes [TT] . 
induced by top and bottom quark loops, have the top 
loops dominant at small tan /3, resulting in a decrease of 
the total cross section with tan /3 up to the point where 
the b loops take over and the total cross section increases. 
Finally, the UH production mode is kinematically sup- 
pressed and decreases with 1/ tan 2 /3, while vector boson 
fusion and associate production with gauge bosons is not 
important, contrary to the case for the lightest Higgs bo- 
son. 

The decay A/H — »• r + r~ is the main process for the 
LHC experiments to search for the neutral heavy Higgs 
bosons at the present LHC energy, the dominant decay 
into bb being overwhelmed by the SM multi-jet back- 
ground. As such, the tt mode has so far attracted most 
of the attention in the LHC searches for heavy Higgs 
bosons. At intermediate to large values of tan /3 the tt 
and bb channels saturate the decay widths of the A and 
H. At low tan /3 the decay pattern of the heavier MSSM 
Higgs particles becomes more complicated by the onset 
of several decay modes which compete with tt, in par- 
ticular WW, ZZ, tt and hh. The branching fractions for 
the decays of H and A bosons are shown in Fig. [2] as a 
function of tan/3 for two masses below (300 GeV) and 
above (500 GeV) the tt threshold. 

The main features can be summarised as follows. Be- 
low the tt threshold, the H boson decays into gauge 
bosons H — » WW, ZZ and into pairs of the light 
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FIG. 2: Decay branching fractions for H (left) and A (right) 
bosons of mass 300 GeV (above) and 500 GeV (below) as a 
function of tan/3, in absence of decay channels into SUSY 
particles. 



Higgs boson, hh have substantial rates. In the interval 
2Mh < Mjj < 2m t , this interesting channel, H — >• hh, 
becomes the dominant decay mode for tan /3 ~ 3. Sim- 
ilarly the pseudoscalar A boson decays into hZ with a 
significant rate, above threshold and at low tan/3. For 
heavier H/A masses, the top decay channels, H^A^r ti, 
becomes dominant below tan f3 ~ 5 — 10. 

Figure [3] shows the regions in the [Ma — tan jS\ plane 
of the pMSSM parameter space, where the branching 
fractions of H -+ ZZ/WW, H/A ti and H -+ hh 
are larger than 1%, 5% and 10%. As can be seen, 
BR(H/A -> ZZ/WW, hh) can be large in the small to 
intermediate tan /3 and Ma region. Above the threshold, 
the BR(H — > ti) is large for tan j3 < 20 independently of 
the Ma value. 

The exclusion limits in the [Ma — tan 0\ plane may be 
modified when some light SUSY particles are present in 
the spectrum. In particular, three scenarios can affect 
the [Ma — tan /3] bounds. First, light SUSY particles, 
with masses < ^M H / A mav induce SUSY decays of the 
H/A states thus reducing the H/A — >• r + r~ branching 
fraction. For Ma < 1 TeV, these SUSY particles can 
be light neutralinos or charginos, indicated collectively 
with x m the following, and light sleptons, in particular 
staus, f , in the case of the CP-even H boson while for 
the pseudoscalar A, only decays A — » f\T2 are allowed. 

The pMSSM, due to the uncorrelated mass values of 
the SUSY particles afforded by its 19 free parameters, of- 
fers a convenient framework for this study, in particular 
by revealing scenarios where decays into SUSY particles 
may be important. The scans used for this study with 
the constraints and the relevant ranges for the variation 
of the pMSSM parameters have been already presented in 



Ref. [121 [T3] . In this analysis, we start from a large sam- 
ple of 2xl0 8 generated pMSSM points and select those 
fulfilling the constraints from LEP data, flavour physics, 
dark matter and g, q direct searches at LHC as discussed 
in Ref. [12]. In particular, we apply the constraints de- 
rived from the rare decay B s —± and direct dark mat- 
ter searches, also providing us with constraints to the 
[Ma — tan j3] parameter space JT4J [15] , and we impose 
that one of the neutral Higgs bosons has a mass in the 
range 121.5 to 129.9 GeV to be consistent with the results 
of the SM Higgs searches, as discussed below. 

The LHCb experiment has recently announced the first 
evidence for the B s — >> /i + \i~ decay and measured its 
branching fraction to be in agreement with the SM ex- 
pectation [16 . This branching ratio is sensitive to the 
Higgs sector, in particular to Ma and tan /?, proportional 
to rsj tan 6 /3/M^ in the large tan/3 limit. Complemen- 
tary information is also obtained by dark matter direct 
detection experiments, in particular the latest XENON- 
100 limits [17], probing the scattering of neutralino with 
matter, which can be mediated by scalar particles. 

The tools used to perform the scans and the analy- 
sis have been presented in Ref. [I4j [18] . Most relevant to 
this study are the calculations of the Higgs decay branch- 
ing fractions and production cross sections. The first are 
computed using the latest version of HDECAY (5.10) [19] . 
The cross section for gg — ^ H/A process is computed 
at NNLO with HIGLU 3.1 [20] [21], that for bb -> H/A 
at NNLO with bbh@nnlo [22] and that for pp bbH 
at LO with HQQ [23 . In addition, we compare the re- 
sults for gg and bb — ^ H/A from these programs to those 
from SusHi [24] and found an agreement within 10-15%. 
The Higgs and superparticle spectra are calculated with 
Softsusy 3.2.3 [25 and Superlso Relic v3.2 [261127] 
computing the dark matter relic density and flavour con- 
straints and providing the central control program inter- 
faced to the other codes. 



B. SUSY Effects in H and A Decays 

There are regions of the MSSM parameter space where 
the rr channel is suppressed and the limits derived in this 
channel are correspondingly relaxed. These may be due 
to direct decays of H/A to SUSY particles or to loop 
corrections to the H/Abb vertices, affecting the H/A —> 
tt branching fraction. 

We consider first the decays of heavy neutral Higgs 
bosons into pairs of SUSY particles. The heavy Higgs 
bosons couple to charginos and neutralinos, primarily to 
identical particles for the mixed gaugino/higgsino states, 
and to different particles in case of pure gaugino or hig- 
gsino states. If the decay to charginos is allowed, it domi- 
nates over the decays to neutralinos. Heavy neutral Higgs 
bosons also couple to scalar fermions. However, decays 
to scalar fermions of the first two generations are sup- 
pressed and only significant at low tan /?, where they are 
sub-dominant. For scalar fermions of the third genera- 
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FIG. 3: Branching fraction for H — »• ZZ (left), —> tt (centre) and — ► hh (right) for the selected pMSSM points in the [Ma —tan 0] 
plane. 



tion the decay rates can be much larger, but they are 
suppressed at large tan j3 for scalar top quarks, while 
they are enhanced for the scalar taus and scalar bot- 
toms. Since the lightest scalar tau, f\ is often the NLSP 
at large tan /3, decays to staus are usually the dominant 
channel for decays into scalar fermions. 

Figure [4] shows the decay branching fraction of H into 
any pair of SUSY particles calculated for the accepted 
pMSSM points for which at least one of these decay chan- 
nels is kinematically allowed. In approximately 25% of 
these cases the branching fraction into SUSY particles is 
larger than 0.10. 




BR(H -» SUSY) 



FIG. 4: Decay branching fractions for H into SUSY particles. 
In the left panel branching fraction for inclusive decays into 
any pairs of SUSY particles for the accepted pMSSM points 
where at least one of these decay channels is kinematically 
allowed. In the right panel, branching fraction into pairs of 
charginos and neutralinos, XX as a function of //tan/?. 

The yield in the XX channels, representing the sum of 
all the kinematically accessible chargino and neutralino 
pairs, depends on the mass parameter M 2 and the Hig- 
gsino mass mixing parameter /i. Figure [5] shows the 
<j x BR product in the [pi — M2] parameter plane to high- 
light the enhancement of this class of decays along the 
small M2 or fi regions. 

The rates of decays into SUSY particles depend mostly 
on the difference between the masses of heavy boson and 
those of the SUSY particles. As the scale of the mass of 
the H and A bosons probed at the LHC increases, de- 
cays into SUSY particles become more likely and have to 
be carefully considered. The relevant mass patterns are 
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FIG. 5: Product of production cross section and decay 
branching fraction for H — > XX a ^ 14 TeV in the [fi — M2] 
parameter plane. The dots in the light colour show all the 
selected pMSSM points and those in darker shades of colour 
the points having a x BR larger than 1, 10 and 100 fb. 



extensively probed in our pMSSM scans. The increase in 
the branching fractions of any of these SUSY channels is 
correlated to the decrease of that for the tt mode, which 
can be suppressed by a factor of two, or more, compared 
to its average value at large Ma and tan j3 values. 

Finally loop corrections to the Hbb and Abb vertices, 
known as A5 corrections [28], modify both the H/A pro- 
duction rates and their decay widths. In the decoupling 
limit, the H/A coupling to bb is modified by a factor 
(1 + A 6 ) _1 , where 



2as 



/iM 3 



tan /3 



max(M|,m~ , m~ ) 



16tt 2 max(// 2 ,ra~ , m~ 



(i) 



Full one loop corrections to the WW and ZZ decays 
have also been computed [29j [30] . We observe that the 
BR(H/A — » tt) is reduced as a result of the enhancement 
of BR(H/A — ^ bb) due to these corrections for SUSY pa- 
rameters yielding a large A5 of negative sign (see Fig- 
ure |6|. Such a large, negative A5 term has also implica- 
tions on the decay branching fractions of the lightest h 
boson. In the region where the H — >> tt decay rate is re- 
duced, the branching fraction BR(/i — » bb) is also reduced 
and those for the other modes correspondingly increased. 
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These patterns might be tested through more precise de- 
terminations of the signal strengths of the lightest Higgs 
decays. In view of these effects, redundancy obtained 
through search in multiple channels sensitive in the same 
regions of the [Ma — tan /3] parameter space appears to 
be essential. 




FIG. 6: Correlation of the BR(H -> rr) with BR(H -> bb) 
(left) and with A& (right) for accepted pMSSM points. The 
correlated suppression of both the rr and bb branching frac- 
tions are due to additional decays into SUSY particles, while 
the decrease of rr with the increase of bb is due to A& effect. 



C. Constraints from the ft Mass and Decay Rates 

Assuming that the observed ~126 GeV state is the 
lightest Higgs boson of the MSSM, ft, its mass Mh de- 
pends on several SUSY parameters, in particular Ma, 
tan/3 and the SUSY scale, M s = y/m^m^. The LEP-2 




M A (GeV) 

FIG. 7: The [Ma — tan/3] parameter space compatible with 
121.5 <M h < 129.9 GeV for different SUSY scales M s . Dis- 
tribution of the accepted pMSSM points in the [Ma — tan /3] 
compatible with the Mh mass interval for 0.5< Ms <3.5 TeV 
(black dots) and 5< M s <20 TeV (light grey dots). 

limit [31], Mh > 114.5 GeV has been long used to define 
a constraint in the [Ma — tan j3) plane, corresponding to 
tan/3 > 2.4 for M s = 1 TeV and M top = 172.9 GeV [32 
in the so-called M™ ax scenario [33]. Now, each value of 
Mh defines a constraint in the [Ma~ tan /3] which depends 
on Ms- A larger value of Ms corresponds to a weaker 
constraint on tan /3. Therefore, it is possible to set a large 
enough Ms scale which recovers the low tan /3 solutions 
of the MSSM, even by applying the LHC constraints for 



Mh, which are now significantly stronger than the LEP- 
2 limit. This is illustrated in Figure [7] which shows the 
values of tan /3 vs. Ma for our pMSSM scans which are 
compatible with 121. 5< M h <129.9 GeV, for two inter- 
vals of values of the SUSY scale Ms- Large enough Ms 
values rescue the MSSM scenarios at low values of tan/3, 
provided we accept a high fine tuning parameter from the 
large scale of Ms. This observation motivates the special 
attention we have chosen to devote to low tan (3 scenarios 
in this study. 
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FIG. 8: Scaling of the ft branching fractions into bb, rr and 
ZZ normalised to their SM value as a function of tan /3 for 
Ma = 300 GeV (left panel) and as a function of Ma for tan /3 
= 10 (right panel). 



A second set of indirect constraints is derived by the 
measured ft decay rates. For large Ma values, the cou- 
plings of the ft boson can be expanded in powers of 
Mz/Ma to obtain the following tree- level result [TO] : 



Ma>M z , 
9hVV — > 1 



M%_ 



sin 



2 ^ tan/3»l ^ 



2M| 



M\ tan 2 /? 



(2) 



For Ma ^> Mz, gnvv reaches the SM value, more quickly 
if tan /3 is large. The ft couplings to up- and down- type 
fermions scale as [TO] : 



9hz 



Qhdd 



M A ~>M Z 



M A ^M Z 



Mf Sin 4/3 tan/3>l 



2M 2 A tan/3 



2M| 
M 2 A tan 2 /3 



(3) 



Mf 
2M\ 



sin 4/3 tan /? 



tan/3>l 



M (4 ) 



The couplings of the ft boson approach those of the 
SM Higgs boson for Ma ^> Mz and these limits are 
reached at lower values of Ma for large tan/3 (see Fig- 
ure |8|. In practice, the ratio of branching fractions 
Rxx = BR(ft -> XX)/BR(# S M XX) or the sig- 
nal strengths nxx — / &(H$m) x Rxx, where a is 
the relevant production cross section, can be used to set 
constraints on the value of Ma- The recent approximate 
N 3 LO calculation of the Higgs production cross section 
resulting in a 17% correction also needs to be taken into 
account [34] . The latest set of LHC results already allows 
us to evaluate some non- trivial constraints, as discussed 
in the next section. 
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III. CONSTRAINTS IN THE Ma — tan PLANE 

The LHC searches have gathered a significant corpus 
of results, which can be used to place some important 
constraints on the H and A bosons in a variety of chan- 
nels. These results also allow us to study the expected 
sensitivity of data to be taken at 13 and 14 TeV from 
2015. In the next two sections we discuss the current 
constraints and in the following we present the extrap- 
olation to 14 TeV. There are important decay channels, 
such as H/A — >• hh and A — >> hZ, for which no analy- 
sis has been performed yet on the LHC data. We char- 
acterise the kinematics and reconstruction strategy for 
these processes using parametrised simulation at the end 
of this section. 



A. Present Constraints (7 and 8 TeV) 

1. Indirect Constraints from the light Higgs signal 

The ATLAS and CMS collaboration have recently 
updated their determination of the mass and signal 
strengths of the Higgs-like particle. In particular, re- 
sults for the 77 [33 [36], [33 [38] and WW [39] [40] 
channels have been reported by both collaborations for 
the full 8 TeV data set, corresponding to integrated lumi- 
nosities of up to 25 fb _1 . In addition, CMS has updated 
the search in the rr channel at low mass [41] . Here we use 
the weighted averages for the mass and signal strengths 
in these preliminary results, as summarised in Table [l| 
For the important 77 channel we average the prelimi- 
nary results of the multi-variate and cut-based analyses 
of CMS accounting for the quoted correlation [48 . The 
results of the two collaborations are only marginally con- 
sistent and we therefore rescale the error of the com- 
bined result according to the prescriptions of the Parti- 
cle Data Group [42]. We use these new inputs and per- 
form an analysis of the regions of MSSM parameter space 
favoured by these data. The analysis follows the strategy 
discussed in [13]. We define the 90% C.L. region corre- 
sponding to observables given in Table [I] by constructing 
the corresponding x 2 probability. We account for the the- 
ory uncertainties on the MSSM h mass, ±1.5 GeV, and 
the Higgs production rates, ±20%. No signal evidence 
has been reported for the 66, where we also include the 
combined estimate on ^55 obtained by CDF and DO at 
the Tevatron [45 , and the r + r~ channels. For these we 
add the contribution to the total % 2 only when the re- 
spective ji value is outside the ±1.5 a interval from the 
measured central value. Compared to the results avail- 
able at the end of 2012, we register a marked realignment 
of the average values for the jixx signal strengths from 
the ATLAS and CMS results around the SM values. This 
has important consequences on the constraints derived. 
In particular, the Ma bound derived from the new data 
is about 100 GeV lower compared to that obtained on 
the first preliminary results on part of the 8 TeV data re- 



leased at the end of 2012, without the CMS re- analysis of 
the 77 channel [13] . This clearly shows that it is difficult 
to predict the sensitivity achievable in future for these 
indirect limits, which depends not only on the accuracy 
of the inputs but also rather critically on the measured 
values. 



Parameter 


Value 


Experiment 


M h (GeV) 
Hzz 

flWW 


125.7±0.4 
1.20±0.30 
1.10±0.22 
0.77±0.21 


ATLAS 46 +CMS[38 
ATLAS[35]+CMS[36] 
ATLAS[37]+CMS[38] 
ATLAS 39 +CMS[40 


fl TT 


1.12±0.45 
1.01±0.36 


ATLAS 43 +CMS[44 +(CDF+D0)[45 
ATLAS [47] +CMS[4T] 



TABLE I: Input values for the average values of the h mass 
and signal strengths used for this study with their statistical 
accuracies. Systematic uncertainties are discussed in the text. 




150 -3 -2-10123 
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FIG. 9: Fractions of pMSSM points compatible at 90% C.L. 
with the constraints of Table [I] in the [Ma — tan f3] (left) and 
\p,-M 2 ] (right). 

We consider points compatible at 90% C.L. with these 
inputs accounting for theory uncertainties. We observe 
that these account for 76% of the accepted pMSSM 
points, up from the 30% obtained in the same analy- 
sis performed on the earlier data. For all these points 
the -126 GeV state observed by ATLAS and CMS is the 
lightest Higgs, h. Therefore we confirm the results from 
our previous analysis where we did not find any pMSSM 
solution compatible with the LHC Higgs results where 
the 126 GeV particle is either the H or the A boson. 
This result provides an answer to the question of [49J. 
Figure [9] shows this fraction as a function of [Ma — tan /?] 
and [/i — M2]. From the [M^ — tan/3] distribution of these 
points we define the region containing 99% of the points 
compatible at 90% C.L. with the LHC Higgs results. This 
region defines an indirect lower bound on Ma-, which will 
be compared to the direct exclusion from H/A searches 
in the next section. 



2. Direct Constraints from MSSM Higgs Searches 

Searches for the H/A —> r + r~ process have been con- 
ducted by the ATLAS with 4.7 fb" 1 at 7 TeV [50 and 
CMS with 4.8+12.2 fb" 1 at 7 and 8 TeV [EE]. The CMS 
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M(A) (GeV) M(A) (GeV) M(A) (GeV) 

FIG. 10: Product of production cross section and decay branching fraction for H — »• ZZ (left), H — > WW (upper centre) and 
H —> ti (right) at 8 TeV in the [Ma — tan 0] parameter plane. The dots in the light colour show all the selected pMSSM points 
and those in darker shades of colour the points having a x BR larger than 1, 10 and 100 fb. The lines superimposed on the left 
panel show the expected (dashed) and observed (continuous) 95% C.L. upper limits obtained in the H/A — > tt search of |51j . 
Entries below threshold in the tt channel are due to off-shell decays. 



sensitivity corresponds to an expected upper limit on the 
product of production cross section and decay branching 
fraction of -80 fb at 300 GeV and 20 fb at 500 GeV. 
In this study, we impose the expected CMS 95% C.L. 
limit on the product of production cross section and de- 
cay branching fraction, which is weaker than the observed 
limit, on our pMSSM points. 

The production and decay pattern of the heavy MSSM 
neutral Higgs bosons crucially depend on the value of 
tan /3, as discussed above. The LHC data at 7 and 8 TeV, 
probe relatively large values, tan/3 > 5-10. For these 
values, their couplings to b quarks and r leptons, propor- 
tional to tan /3, are strongly enhanced, and those to top 
quarks and massive bosons, proportional to ~ 1/ tan/3, 
are suppressed. Therefore the rr channel is presently 
the single most constraining decay mode. It defines a 
region of the [Ma — tan /3] parameter space which is 
probed also by the B s — » jiji rare decay and by dark 
matter direct detection experiments p~4j [T5] , but the 
H/A —> tt LHC searches at 7 and 8 TeV set the tightest 
constraints. In addition to it, preliminary results have 
been reported for the first search for H/A — >• bb in as- 
sociate production with b jets bbH / A — ^ bbbb based on 
the 7 TeV CMS data [52 , which has sensitivity at large 
values of tan j3 with an expected upper limit of 8 pb on 
a x BR at 300 GeV. The analyses of the decays of the SM 
Higgs Hsm — > ZZ have set constraints on the product 
of production cross section and decay branching fraction 
<J (99 ~ > Hsm) x RR(Hsm ZZ) for Higgs masses up 
to 1 TeV [53l [54] . These can now be used to constrain 
the decay of the heavy SUSY H/A —> ZZ, with upper 
limits of -1.9 and 1.4 fb at 200 and 300 GeV, respec- 
tively. Finally, the decay H/A — >• tt can be constrained 
through the cross section bounds obtained for the pro- 
duction of a narrow resonance decaying into top quark 
pairs, interpreted in the original studies in the context 
of the searches for the production of a lepto-phobic Z' 
gauge boson, KK resonances or other exotic narrow reso- 
nances. Results have been reported by both ATLAS [55] 
and CMS [56] for the 7 TeV data with cross section upper 
limits of order of 3 pb and 0.8 pb at resonance masses of 



-500 GeV and 800 GeV, respectively. 

First, we study the value of the product of produc- 
tion cross section and decay branching fraction in sev- 
eral channels by scanning over the pMSSM parameters. 
Figure 10 shows the regions of the [Ma — tan /3] param- 
eter space where the product a x BR exceeds 1, 10 and 
100 fb at 8 TeV for the gg H/A and bb H/A pro- 
duction processes and the H/A -» ZZ, H/A ->> WW 
and H/A — >• tt decays. Finally, we combine the con- 
straints derived in the various channels. We take the 
expected upper limits on the products a x BR in the 
various channels for both (a) the present status of the 
results and (b) their extrapolation to the full 8 TeV data 
set of 25 fb -1 . When limits are only available for the 
7 TeV data set, we compute the expected limit at 8 TeV 
by taking the ratio of production cross sections at the 
two energies, as a function of the H/A mass, into ac- 
count. For each channel, we consider the contours in the 
[Ma — tan /3] plane where more than 95% of the selected 
pMSSM points are excluded by these constraints. Along- 
side the rr channel, the ZZ and bbbb channels also offer 
sensitivity on the 7 and 8 TeV data. For the ZZ chan- 
nel we use the upper limits on a x BR from [54]. These 
limits define an excluded region which connects with the 
rr constraint at low masses and extends up to Ma — 
550 GeV for tan/3 =3-4. The bbbb channels is based on 
the preliminary result on the 7 TeV CMS data [52] ex- 
trapolated to the full 8 TeV data set. We include also 
the constraint derived from the signal strengths, \i, ob- 
tained in the ATLAS and CMS SM Higgs analyses for 
the 77, WW, ZZ channels and the limits for bb and rr, 
as discussed above, interpreting the observed particle as 
the SUSY lightest Higgs, h. Results are summarised in 
Figure 11 The combination of the H/A —> tt channel 
and the mass and \i values for the lightest h boson ex- 
clude the region with Ma > 320 GeV for all values of 
tan/3. For the current results, the \i values defines this 
lower bound in the region of tan/3 = 2-15, where the di- 
rect search sensitivity is weaker. The sensitivity of the 
direct H/A searches should approach this bound down 
to tan /3 ~ 10, once the full 2012 data is analysed. The 



8 




1500 

M(A) (GeV) 




1500 

M(A) (GeV) 

FIG. 11: Combination of the expected constraints on the 
[Ma — tan /3] parameter plane from the rr and ZZ chan- 
nels for (a) the current results (upper panel) and (b) their 
extrapolation to the full 8 TeV data set (lower panel). The 
colour scale gives the fraction of pMSSM points excluded at 
each Ma and tan ft value. The contours show the limits cor- 
responding to 95% or more of the points excluded. The 90% 
C.L. constraint from the Higgs signal strengths is also shown. 
The expected and observed upper limits on tan /3 obtained in 
the MSSM M™ ax scenario from the rr channel search of [51] 
are indicated by the grey dotted and continuous lines, respec- 
tively, on the upper plot. The grey region has no accepted 
pMSSM points after the B s — > direct DM searches and 
Mh constraints. 



ZZ channel, and to a lesser extent the WW, should close 
the low Ma corner from tan (3 ~ 2 up to the rr limit for 
M A < 230 GeV with the full 8 TeV data. The upper 
limits from the tt channel, for which only results at 7 TeV 
have been reported, are still below the expected values 
for H/A production in the MSSM, even by extrapolating 
them to the full 8 TeV data set. Instead, this channel 
will become essential at 13 and 14 TeV. The combina- 
tion of these constraints from the Higgs sector provide 
limits on Ma and tan /3, which are significantly tighter 
compared to those derived from flavour physics, such as 
the BR(£? S — » jiji) for which the first measurement has 
recently been reported by LHCb [16] (see Figure 11). 



B. Perspectives at 14 TeV 

The increase of the production cross sections moving 
from 7 to 14 TeV is a factor of 4.5 to 9 for gg -> H/A and 
5 to 12 for bb H> H/A in the mass range 300 to 800 GeV. 
! shows the regions of the [Ma - 



Figure 
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-tan/3] parameter 
space where the product a x BR exceeds 1, 10 and 100 fb 
for the gg — > H/A and bb —> H/A production processes 
and the H/A -> ZZ, H/A WW. H/A -> hh, A -> 
hZ, H/A -» tt and the inclusive decays H/A -> SUSY 
particles. At the high mass end the product a x BR 
of ~10 fb, corresponding to the current sensitivity at 
800 GeV in the rr channel, is obtained beyond Ma = 
1 TeV. At 13 and 14 TeV the sensitivity extends to mass 
values above the hh, hZ and the tt production thresholds 
at small to intermediate values of tan /3, which make these 
channels relevant to the LHC searches. In this region 
the rr channel alone cannot ensure the coverage of the 
[Ma — tan fS\ plane and these additional channels need 
to be included. The ZZ channel provides redundancy 
while the tt decay is most important, in particular at 
large Ma and low tan f3 values. The WW channel has 
more limited interest, since its sensitivity is lower than 
ZZ. The combination of the rr, ZZ and tt modes covers 
the [Ma — tan f3] parameter plane up to Ma — 700 GeV 



for any value of tan /?, as shown in Figure 13 



1. Characterisation of hZ and hh channels 

The decays H —> hh and A — >• hZ are important in 
providing redundancy at low values of tan f3 and interme- 
diate Ma masses. They also result in rather distinctive 
bbbb, bbrr and bbit (£ = e, ji) final states, which should 
be investigated in the high energy LHC runs. Since these 
modes have not yet been searched for in the LHC data, 
we characterise here their decay kinematics and study 
the reconstruction strategies using a simple analysis for 
signal events. 

These events are generated using Pythia 8.1 [57] 
at 14 TeV and scaled to an integrated luminosity of 
150 fb -1 . For this study, we have chosen Ma = 400 
and 500 GeV, tan/3 = 5 with branching fractions of 
0.12 for H -> hh and A -> hZ. The detector re- 
sponse simulation is performed using Delphes 3.0 [58] . 
Jets are reconstructed using the anti-kt algorithm [59] 
implemented in Fast Jet [60] . requiring their pseudo- 
rapidity, 77, not to exceed 2.8 and transverse momen- 
tum pt >20 GeV. Electrons and muons are accepted 
for \rj\ <2.4 and pr >20 GeV. 6-jets are accepted at 77 < 
2.5, assuming a tagging efficiency of 75% per jet. In both 
channels, b jets are rather soft, with the transverse energy 
distributions peaking around 50 GeV, thus emphasising 
b tagging at relatively small transverse energies (see Fig- 
ures 14 and 15). Similarly low is the transverse energy 



distribution of leptons from the Z decay in the A chan- 
nel, which has its most probable value just above the pt 
cut applied in this analysis (see Figure 15). 



9 




800 1000 1200 1400 
M(A) (GeV) 



600 800 1000 1200 1400 
M(A) (GeV) 




Selected Points 

a x BR(H -» hh)> 1 fb 

ct x BR(H — » hh)> 10 fb 

axBR(H-»hh)>100fb 



200 400 600 800 1000 1200 1400 
M(A) (GeV) 



Selected Points 

o x BR(H -» hZ)> 1 fb 

oxBR(H-»hZ)>10fb 

oxBR(H-»hZ}>100i 




1200 1400 
M(A) (GeV) 



Selected Points 
ct x BR(H — » tt)> 1 fb 
ct x BR(H —* tt)> 1 f b 
a x BR(H -» tt)>100 fb 




1200 1400 
M(A) (GeV) 



Selected Points 




200 400 600 800 



1000 1200 1400 
M(A) (GeV) 



FIG. 12: Product of production cross section and decay branching fraction for H — >> ZZ (upper left), H — >• WW (upper centre), 
H — » hh (upper left), A — > hZ (bottom right), H —> ti (bottom centre) and H — > SUSY particles (bottom left ), at 14 TeV in 
the [Ma — tan 0] parameter plane. The colour coding is given in the legend and it is the same as in Figure 10 



H — >• hh — >• bbbb events are reconstructed by requiring 
at least three 6-tagged jets. The pairing of four b jets, or 
three b jets with any of the reconstructed jets, which min- 
imises the mass difference of the two di-jet pairs and their 
difference from the h mass of 126 GeV is selected. The 
di-jet invariant mass distribution is shown in Figure 14 
The invariant mass resolution obtained with the fast sim- 
ulation is comparable to that reported for the H^m bb 
search. The four-jet invariant mass, M bbbb shows a clear 
peak corresponding to the generated H mass as shown in 
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FIG. 13: Combination of the expected constraints on the 
[Ma — tan /3] pa ram eter plane from the rr, ZZ and tt chan- 
nels as in Figure [TT| extrapolated to 150 fb" 1 at 14 TeV. The 
colour scale gives the fraction of pMSSM points excluded at 
each Ma and tan /3 value. The grey region has no accepted 
pMSSM points after the B s — >> fifi, direct DM searches and 
Mh constraints. 



Figure 14 The efficiency of this selection for the signal 
mass region of 300< M bbbb < 500 GeV is ~ 16% at both 
values of Mh- 




FIG. 14: Reconstruction of H -> hh -> bbbb events at 14 TeV 
for Mh — 400 GeV: distribution of the fr-jet transverse energy 
Et (upper right) and energy E (lower left), invariant mass of 
bb pairs (lower left) and bbbb invariant mass (lower right). A 
BR(if — > hh) = 0.12 has been assumed. 

For the Zh — >• tibb we select events with two, oppo- 
sitely charged, electrons or muons with two or more jets, 
of which at least one b tagged. The ££ invariant mass 
is required to be consistent with that of the Z within 
the resolution. If the event contains exactly two 6-tagged 
jets, the invariant mass of the pair is required to be con- 
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FIG. 15: Reconstruction of A Zh ££66 events at 14 TeV 
for Ma = 400 GeV: distribution of the lepton transverse 
energy pr (upper left), 6-jet transverse energy Et (upper 
right), bb (lower left) and bbll (lower right) invariant mass. 
A BR(A -> Zh) = 0.12 has been assumed. 



sistent with 126 GeV within the resolution. If there is 
only one 6-tagged jet, but it has a mass consistent with 
126 GeV, this is also accepted. The final mass is com- 
puted by combining the di-leptons with the di-jet pair 
or the single b jet. The resulting distribution is shown 
in Figure 15 for an integrated luminosity of 150 fb _1 . 
The selection efficiency for the loose signal mass region 
of 300< M bb u < 500 GeV is ~ 25% at both values of M A . 




1500 

M(A) (GeV) 

FIG. 16: Regions on the [Ma — tan/3] parameter plane where 
the H — >• hh and A — > hZ process yield 50 reconstructed 
events for 150 fb _1 at 14 TeV, compared with the coverage 
provided by the combination of rr and ZZ shown in dark 
blue. The grey region has no accepted pMSSM points after 
the B s — > /i/i, direct DM searches and Mh constraints. 

Since we base this preliminary characterisation on the 
reconstruction of signal only event and have not consid- 
ered the backgrounds, we cannot define here exclusion 
contours. Instead, we simply plot the regions of the 
[Ma — tan f3] plane where we register more than 50 recon- 
structed events for 150 fb _1 of data at 14 TeV. The result 
is shown in Figure [16} where the region covered by the 



hh and hZ final states is compared to that of expected 
sensitivity for the combination of the rr, ZZ and tt chan- 
nels, considered above. We notice that the hh channel 
covers the full tan f3 range of interest from threshold up 
to M A ~ 400 GeV and up to 550 GeV at low tan j3 val- 
ues, beyond the ZZ sensitivity. In this important region 
of small to intermediate values of tan/3, the hh and hZ 
channels provide redundancy to the coverage offered by 
the rr and it modes. 



Effect of QCD Uncertainties and SUSY 
Particles 



The limits derived above do not account for the ef- 
fects of theoretical uncertainties, affecting the Higgs pro- 
duction cross section and decay branching fractions, and 
of SUSY contributions. First, the gg — » H/A and 
bb — » H/A cross sections have sizeable QCD uncertainties 
from the factorisation and renormalisation scales, par- 
ton distribution functions (PDFs) and parametric sys- 
tematics from a s and the heavy quark masses. We es- 
timate the parametric systematics on the cross section 
for a s = 0.118 ±0.0012, m h {m h ) = (4.19±0.05) GeV, m t 
= (172.9±1.5) GeV and those from the PDFs by taking 
the largest difference between different sets of functions. 
The latter is the dominant contribution. The combina- 
tion of the uncertainties on the quark masses, PDFs and 
a s leads to an estimated systematic uncertainty on the 
pp -> H/A rate of w ±24% at 8 TeV and « ±20% at 
14 TeV, dominated by the PDFs and scale, and compa- 
rable to those for pp — > Hsm production [61] [62] . 

In order to evaluate their impact on the exclusion con- 
tours in the [M A — tan /?] plane, we repeat our study while 
changing the production cross section by ±25% and com- 
pare the constraints obtained to that corresponding to 
the central values for the production cross sections. Fig- 
ures [TT] and [18] show the fractions of excluded points in 
the [Ma — tan /?] plane and their projections as a function 
of Ma for the fixed value of tan f3 = 15 at 8 and 14 TeV, 
respectively, and includes the effect of the ±<jqcd change 
of the cross sections by the QCD uncertainties. The ef- 
fect is a shift of the excluded Ma mass by ±45 GeV at 
8 TeV and by ±55 GeV at 14 TeV at tan /3=15 and larger 
for higher values of tan j3. 

Then, we observe that, there is a significant smear- 
ing of the curve giving the fraction of excluded pMSSM 
points as a function of Ma, even if the systematics on the 
production cross section are ignored. In fact, the exclu- 
sion curve goes from 10% to 90% of the points excluded 
over a range of Ma values spanning ~90 GeV at 8 TeV 
and ~150 GeV at 14 TeV, as a result of the variation 
of other pMSSM parameters. This range, which is com- 
parable to that corresponding to the QCD uncertainty 
obtained above, is intrinsic to the pMSSM and includes 
contributions such as the loop effect through the A5 term 
discussed in section II. B. 

Finally, we consider quantitatively the region of the 
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FIG. 17: QCD systematics and SUSY particle effects on the 
projected H/A — > rr exclusion with at 8 TeV. Left: fraction 
of pMSSM points excluded for tan ft — 15 as a function of 
Ma (continuous line) and the effect of a change by ±25% 
of the pp — »• H/A production cross section to reflect QCD 
uncertainties (dashed lines). Right: limits from the rr in 
the [Ma — tan fi\ plane obtained by varying the production 
cross section (dashed lines) and requiring less than 0.1% of 
the points around the limit to fail exclusion due to the effect 
of SUSY decays. The grey region has no accepted pMSSM 
points after the B s — ► fifi, direct DM searches and Mu con- 
straints. 
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FIG. 18: QCD systematics and SUSY particle effects on the 
projected H/A — » rr exclusion, as in Figure 17 for 150 fb _1 
at 14 TeV. 



[Ma — tan /3] plane where decays into SUSY channels may 
invalidate the rr limit. The panels on the right in Fig- 
ures f 



17 and [181 nave clashed lines showing the limit of the 



region where this may occur, for 8 and 14 TeV, respec- 
tively. Since the axBK(H/A — >• rr) product increases 
for low Ma and high tan f3 values, there is a region where 
the SUSY decays cannot upset the exclusion obtained in 
this channel, since the H — » SUSY branching fraction is 
< 0.60, as shown in Figure [4| However, the region af- 
fected by the SUSY decays expends much further towards 
lower A masses compared to that describing the effect of 
the QCD and parametric uncertainties on the production 
cross section and also the SUSY loop effects. The width 
of this region, ~150 GeV at 8 TeV for tan /? =15, and the 
occurrence of these points increase with the energy which 
gives access to heavier bosons with decays into pairs of 
SUSY particles kinematically allowed. Moving from 8 to 
14 TeV, the width of the regions doubles and the occur- 
rence of these points increase by a factor of ~ 1.5. The oc- 
currence of the various SUSY effects we have discussed for 
the pMSSM points within the region of the [Ma — tan j3] 
plane excluded in the MSSM M™ ax model but not ex- 
cluded due to their low value of the axBR(H/A — >> rr) 
product has been studied for our scans. We observe that 
H/A decays into XX pairs are responsible for 55% of the 



cases and those into ff pairs for another 10%, while it 
is a A5 term large and negative in sign which suppresses 
the rr rate for the remaining 35% of the points failing 
exclusion. When the H decays into XX pair, the domi- 
nant state is the lightest neutralino Xi i n about 1/4 of 
the cases. In the other cases, we observe an increase in 
the yield of X2 3 ~> V^Xi? which may offer an important 
signature for the LHC searches [63] . 



IV. CONCLUSIONS 

The search for heavy Higgs bosons represents a next 
frontier in the understanding of the Higgs sector after 
the discovery of the Higgs-like state at 126 GeV at the 
LHC and the first results on its decays, spin and parity. 
The combination of the indirect limits from the h signal 
strengths and the direct searches in the rr and ZZ chan- 
nels should impose an exclusion limits in [Ma — tan f)\ 
plane around M A > 320 GeV for the 7+8 TeV data, 
determined by the indirect limit from the rates of the 
observed Higgs boson 

As the mass sensitivity of the LHC searches increases 
with the energy and integrated luminosity, more final 
states than rr become relevant to effectively constrain 
the supersymmetric parameter space, in particular at low 
to moderate values of tan/3. In fact, low values of tan/3 
are still viable, after incorporating the M^ constraint, 
provided high SUSY scales, Ms, are chosen and they 
represent a scenario, rich in decays into tt and ZZ, hh 
and hZ boson pairs, which should be carefully explored 
at the LHC at 13 and 14 TeV. The M s bound will reach 
M A > 800 GeV for any value of tan/3 with 150 frr 1 
of data at 14 TeV, determined by the direct searches 
for heavy Higgs states. The effects of the SUSY particle 
spectrum, other SUSY parameters and the QCD theoret- 
ical uncertainties need to be carefully considered. SUSY 
loops and QCD effects on the Ma bounds are find to 
be quite comparable in size. However, scenarios where 
decays into SUSY particles are important, or even dom- 
inant, exist and these channels need to be accounted for 
in the LHC searches at 13 and 14 TeV. 

The constraints derived by the study of the Higgs sec- 
tor are becoming an essential part of the probe of the 
SUSY parameter space at the LHC and offer an essen- 
tial complement to the searches for strongly interacting 
SUSY particles and gauginos. 
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